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Abstract 

The Euler's equations describe the motion of inviscid fluid. In the 
case of shaUow water, when a perturbative asymtotic expansion of the 
Euler's equations is taken (to a certain order of smallness of the scale 
parameters), relations to certain integrable equations emerge. Some 
recent results concerning the use of integrable equation in modeling 
the motion of shallow water waves are reviewed in this contribution. 
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1 Governing equations for the inviscid fluid mo- 
tion 

The motion of inviscid fluid with a constant density p is described by the 
Euler's equations: 

dv 1 

_ + (vV)v = --VP + g, (1) 
V-v = 0, (2) 

where v(x, y, z, t) is the velocity of the fluid at the point (x, y, z) at the time 
t, P is the pressure in the fluid, g = (0, 0, —g) is the constant Earth's gravity 
acceleration. 

Consider now a motion of a shallow water over a flat bottom, which is 
located at z = 0. We assume that the motion is in the x-direction, and 
that the physical variables do not depend on y. Let h be the mean level of 
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the water and let r]{x,t) describes the shape of the water surface, i.e. the 
deviation from the average level. The pressure is 

P = PA + pg{h-z)+p{x,z,t), (3) 

where Pa is the constant atmospheric pressure, and p is a pressure variable, 
measuring the deviation from the hydrostatic pressure distribution. On the 
surface z = h + r], P = Pa and therefore p = rjpg. Taking v = (n, 0, w) 
we can write the kinematic condition on the surface as [e.g. see (Johnson 
1997)] 

dr] drj 

w = — — h u— on z = n + n. (4) 
ot ox 

Finally, there is no horizontal velocity at the bottom, thus 



w = on z = 0. (5) 
The equations ([I]) - ([5]) give the system 



1 

Ut + UUx + WUz = Px, 

P 

1 

Wt + UWx + WWz = --Pz, 
P 

Ux + Wz= 0, 

w = r]t + ur]x, p = ripg, on z = h + r], 
w = on z = 0. 



(6) 



Let us introduce now dimensionless parameters e = a/h and 6 = h/\, 
where a is the typical amplitude of the wave and A is the typical wavelength 
of the wave. Now we can introduce dimensionless quantities, according to the 
magnitude of the physical quantities, see (Johnson 1997, 2002) for details: 

X — > Ax, z zn, t — > r^ t, rj — > arj, 

Vgri 

u — > e^/ghu, w — > e6y/ ghw, p — > epgh. 

This scaling is due to the observation that both w and p are proportional 
to £ i.e. the wave amplitude, since at undisturbed surface (e = 0) both w = 
and p = 0. The system ([6|) in the new, dimensionless variables is 
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ut + e{uux + wuz) = -Px, 
5'^{wt + e{uwx + WW;,)) = -pz, 

Ux+Wz = 0, 

w = rjt + eurjx, p = 1], on z = 1 + er], 
w = on z = 0. 

(7) 

For the right-running waves one can introduce the so-cahed far field 
quantities, see (Johnson 1997, 2002, 2003a) 



( = ^(x-t), T = e^''^t, w = ^eW, (8) 

and the system ([7]) acquires the form 

eur — + e{uu(^ + Wuz) = —P(, (9) 

e5\eWr -W^ + e^uW^ + WW^)) = -Pz, (10) 

u^ + Wz = 0, (11) 

W = er]r — V( + £:u7]i^, p = r], on z = 1 + er], (12) 

W = on z = 0. (13) 



2 Asymptotic expansion of the variables 

Following the idea of Johnson (2002), we can express the variables u, W, p 
as double-asymptotic expansion (in e and 6) with terms, depending only on 
rj(x, t) and explicitly on z. As a result, a single nonlinear equation for r] will 
be obtained, and thus all variables will be expressed through the solution of 
this equation. 

From (jlOp it is evident that pz = 0{e5'^), and thus in the leading order 
p does not depend on z, i.e. 

p = rj. (14) 

Substitution of (dH) into ([9]) and ([TT]) gives for the leading orders 

u = 7], W = —zrj(^. (15) 

Consider the next terms (first corrections) in the expansion of u and W , 
denoted by E(u), E{W), which possibly contain terms of orders e and (5^. 
Writing u = rj + E{u), W = —zr](^ + E{W), from ([9|) it follows 

E^{u) = eirjr + vvc)^ (16) 
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and from (fTTIl . (fT3]) and (fTBI). 



E{W) = - j E(^{u)dz = -ezi'qr + md^ 

W = -z{r](; + er]r + erjrji^). (17) 



Now the substitution of (jTTj) into ([121) gives the leading order equation 
for 7]: 



Vr = -lm + 0{e,5''). (18) 

From (jlSp and (jl6|) we obtain E{u) = —erf /A, i.e. no (5^ term is present 
and finally, using (fTTj) and 



^^ = ??-|f?^ ^ = -2;(?/c - |wc)- (19) 

Using (fT5]) in (fTOl) we have = —e6'^zr](^i^. This can be integrated due 
to (jl2p and thus we obtain the next order approximation for p: 

1 _ ^2 

p = rj-e6^^—r,^(;. (20) 



We accomplished the first step, i.e. starting from the leading order (fT^ 
P^ . we obtained the first corrections ([^01) and an equation for r/, 
The next step can be performed in a similar fashion and it gives 



« = + + y)'?CC' (21) 



where r]{x, t) satisfies the equation 

= -^VVC + lev\ - + 0{e\6\e6'). (24) 

We observe, that at the end of each step the equation for r/(x, t) contains 
terms of smaller order than those, which appear in the expressions for u and 
W. Since we need an equation, containing terms of order 0{e'^,6'^,e5^), we 
need to perform several intermediate sub-steps, (like p^ . (fTTll ) of the next 
step, which leads to the desired equation 
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-r?r?C + -er] r]c^ 



16 
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1 



Now, we can invert ([^T|) by specifying u at a specific depth, 2;o 
(0 < ^ 1): defining n = t, zq)) we obtain 



(25) 



(26) 



where 



.1 4 



-l<A<i. 
6 - - 3 



(27) 



Note that in (|26p there is no term of order e^, xhe substitution of 
in ([25|) yields: 
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+0(e^(^^e252^e5^). (28) 



Next, we go back to the original variables, introducing T = y/et, X = ^/ex, 
see ([8]), keeping only the scaling with e : 



T = -r, X = -T + C, 

e e 



(29) 



or d(; = dx, edr = dT + dx- Thus, ([28]) yields 



ZD I 



r/29 , 5_ 

—+QXjuxuxx + -^uuxxx 

+0{e^,e6\e^5^,e^6^). (30) 



Further, we add formally {sd'^ /lUxxT — ^S"^ fJ-uxxr) to the left-hand side of 
()30p . where is an arbitrary real parameter. In the first term we substitute 



iiT = —Ux — ^£uux, according to ([30]) : 



u - -e5 fj,uxx^ ^ 



-e 



3 _ ^2.1 
-Ux - T,^uux + ed [-fi - -juxxx 
2 2 b 

29 , 9 \ _ /5 3 

— + 6A - -fj-juxuxx ^ \q~ 2^ juuxxx 

+0(e^e5^e=^(5^e25^). (31) 
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We observe that (I30p . (I3ip do not contain terms of orders e and e^. Thus, the 
set-up from (Johnson 2002) naturahy leads to the conclusion, that equations, 
containing nonlinearities as those, appearing in the equations of Camassa 
&: Holm (1993), Fokas & Fuchssteiner (1981) [called also CH from now on] 
and Degasperis & Procesi (1999), Degasperis et al. (2002) [DP for short], 
are generalizations of the Korteweg-de Vries equation, containing the next 
order term (e^iJ^) in the expansion with respect to the small parameters e, 
5. 

3 Integrable nonlinear equations 

In this section we start from a known integrable equation and we try to 
write it in a form, in which it matches (|3ip or (j30p . For another approach 
for matching between water waves equations and integrable equations see 
(Dullin et al. 2003, 2004). The CH and DP equations can be written as 

{U - Uxx)t = ioUx -{b+ l)UUx + bUxUxx + UUxxx, (32) 

where U = U{X,T), oj is an arbitrary constant, 6 = 2 for CH and 6 = 3 
for DP. There is no other choice of the constant coefficients in front of the 
nonlinear terms, leading to integrable equations, see (Ivanov 2005a). Let us 
change the variables in (j32p according to 

X^X-vT, T^T, U^U + C, (33) 
where v and C are arbitrary constants. Then (I32p acquires the form 

{U - Uxx)t = [o; - (6 + l)C + v\Ux + (C - v)Uxxx 

-{b+l)UUx + bUxUxx + UUxxx. (34) 

It is now clear that via the transforms (j33p one can achieve arbitrary coeffi- 
cients for the linear terms Ux and Uxxx- Let us now consider the following 
scaling of the variables: 

X -X, T^PT, U ^ jU. (35) 
a 

Then ()34p can be written as 

{U - a^Uxxh = riUx + r2Uxxx - aPlib + l)UUx 

WhibUxUxx + UUxxx), (36) 
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where ri^2 are arbitrary constants. In order to match ([36]) to ()3ip up to the 
given order, we need to make the following identifications: 

29 , 9 ,/5 3 \ 

which are compatible iff 



5(6+1) ,,30-9f, 

^ 9b 72b ^ ^ 



Thus, (|38|) and (j27|) show that (j36|) describes water waves at depth 



/116-10 , , 

^0 = (39) 

i.e. CH (6 = 2) corresponds to zq = ~ 0.71 and DP (6 = 3) corresponds 
to zq = -v/i| ~ 0.80. The scaling coefficients from (j37p are 



and, apparently only the product /?7 is determined, i.e. there is additional 
freedom in the choice of P and 7, one can take, for simplicity, just 7 = 1 
and then, finally. 



0,= [le'/H, f3 = ^4l re'/'<5"\ 7=1- (41) 

V 2 ' ^ ^{b + 1) ' ' ^ ' 

Another equation, which passes the integrability check developed in 
(Mikhailov & Novikov 2002; Sanders & Jing Ping Wang 1998; Olver & 
Jing Ping Wang 2000) and is presumably integrable (although we do not 
have a proof of this fact - the test provides only a necessary condition for 
integrability) is 

{U — Uxx + Uxxxx)t = Tif/x — V2UXXX + ^iUxxxxx 

-UUx + UxUxx -UxUxxxx, (42) 

where ri^2 are arbitrary constants. This equation contains nonlinearities, 
similar to those, appearing in the nonintegrable equations studied by Holm 
& Hone (2003). 
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The scaling (l35|l gives 

{U -a^Uxx + cx^Uix)T = a[iTiUx-a^f3V2Uxxx + a'^PT2U5x 

- a^jUUx + a^p-fUxUxx - a^p-fUxU4x, 

(43) 

The matching between (I43p and (I3ip leads to the following identifications: 



= I (44) 
In a similar way, from (I43p we find (assuming again 7 = 1) 



A = -g, ,„ = J11«0.96; (45) 



The terms with fourth and fifth derivative in (14311 are of orders 



a' = ^e'6\ a'p=^e'6\ a'(3^ = ^e'6\ (47) 

and, therefore are small, in comparison to the other terms. 
Another set of integrable equations is of the type 



Ut + Uxxxxx + 2(66 + l)UxUxx + 4(6 + l)UUxxx + 20bU^Ux = 0, 

(48) 

where one can recover the Caudrey-Dodd-Gibbon equation, (Caudrey et al. 
1976) for 6 = 1/4, the Sawada-Kotera equation, (Sawada & Kotera 1974) 
for 6 = 3/2 and the Kaup-Kuperschmidt equation, (Kaup 1980) for 6 = 4. 

It is a natural question to ask, if (I48p can match ()30p . Applying the 
transformation ()33p to (I48p . we can write it in the form 

Ut = TUx - 4(6 + l)CUxxx - U^x - mCUUx 

-2(66 + l)UxUxx - 4(6 + l)UUxxx - 20bU^Ux, (49) 
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where T = u + 206C^ can apparently be arranged to be an arbitrary constant 
with the help of the free parameter v. The scaling ([35]). applied to ([M]) gives: 



Ut = a/3TUx - 4(6 + l)Ca^pUxxx - ct'U^x - mCa/S-fUUx 

-a^(3^[2{6b + l)UxUxx + 4(6 + l)UUxxx] - 20ba(3j'^U^Ux- (50) 
Apparently we need to make the following identifications: 

4(6 + l)Ca^l3 = -e6^, AObCapj = -e, 4(6 + l)a^f3j = —e^6'^, (51) 
6 2 12 

giving 

a3 = ^r-^, 7 = (52) 

2006C2' '2 ^ ' 

The order of the term U'^Ux is 20ba(3j'^ = it is not small in com- 
parison to the other terms, and therefore cannot be neglected. Thus, there is 
no direct match between (1480 and (I30p . however, there is more complicated 
transformation, given in (Fokas & Liu 1996) [based on the Kodama trans- 
form, (Kodama 1985)] providing the link between the water-wave equations 
and the integrable systems (HH]) . 



4 Water waves moving over a shear flow 

So far we have only considered waves in the absence of shear. Now let us 
notice that there is an exact solution of the governing equations ^ of the 
form u = U{z), < z < h, w = 0, p = 0, r] = 0. This solution is nothing, 
but an arbitrary underlying 'shear' flow. Waves of small amplitude (of order 
e) propagating over this underlying flow are studied by many authors and 
here we will partially follow Johnson (2003a) and Burns (1953). The scaling 
for such solution is clearly 

u \/'gh(^{z) + eu^ , 

and the scaling for the other variables is as before. Thus, from ([6]) instead 
of d?!) in this case we have 

Ut + Uux + wU' + e{uux + wuz) = -Px, 
5^{wt + Uwx + e{uwx + wwz)) = -pz, 

Ux + Wz = 0, 

w = r]t + {U + eu)r]x , p = ri, on z = 1 + er/, 
u; = on z = 0. 

(53) 



9 



The prime denotes derivative with respect to z. 

In what follows we need the propagation speed c of the waves in the 
linear approximation, i.e. in the case when in ()53p it is taken e = (5 = 0. 
This velocity is now not independent on U . Since pz = 0, in the linear 
approximation p = rj. Let us introduce a stream function ■0, such that 
u = ipz and w = —ipx- In the case of linear waves we can assume that 
-0 = (j){z)e^''^^~^\ 7] = rjQe'^^^^~^^\ where A; is a wave number and 770 is a 
constant. From (I53p with e = 5 = we now easily find a relation between 
tl (z) and 4>{z): 

(l)'{u -c)-u'(t> + m = ^, m = -{u{i)-c)m, <A(o) = o. (54) 

The first equation in (|54p can be written as 

dzjj-c~ {U-cf 

and can be integrated directly. Imposing the boundary conditions from (I54p 
we finally obtain the following relation for the speed of propagation c (the 
so-called Burns condition): 

dz 

/ ^ = 1. 55 

Jo [U{z)-cf 

For a nondecreasing function U{z), such that U{0) < U{z) < U{1) there 
are always two solutions: c > U{1) and c < C/(0). In the absence of flow, 
[7 = these two solutions are simply c = ±1. The presentation in the 
previous sections corresponds to the choice c = 1. 

Again, we introduce the far field variables, cf. ([8]) 

^ = ^(x-ct), r = e^/^t, w = V£W, (56) 

and the system acquires the form 

eur + {U - c)uc_ + Wil' + e(nn^ + Wuz) = -pQ, 
e5\eWr + {U- c)W^ + e{uW^ + WW^)) = -Pz, 

U(; + Wz = 0, 

W = erjr + ([/ — c)rji^ + eur](^, p = r], on z = 1 + erj, 
W = on z = 0. 

(57) 
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Now let us concentrate to the simplest nontrivial case: a linear shear, 
U{z) = Az, where A is a constant. We choose vl > 0, so that the under- 
lying flow is propagating in the positive direction of the x-coordinate. The 
condition (j55p gives the following expression for c: 

c=^(^A±V^ + A^y (58) 

If there is no shear (A = 0), then c = ±1. 

Of course, a parabolic distribution U{z) = U{l){2z — z'^) would be more 
realistic, but then the solution of (j55p is not so simple, cf. (Burns 1953). 

The solution of the system ()57p can be obtained as a series in e and 5 
following the method, explained in the previous section. Here we present 
the final result, obtained in Johnson (2003a). The equation for r] is 

C^ + C^ + l C(c4 + 4c2 + l) 2 1 e2 

= - c(c2 + 1) + 2(c2 + l)3 - H^)^ 

~ 3c(/+l)=^ ^'^'^^^''' + 4c^ + lie' + Q)VCVCC + (c" + + 3)w^cc] 

+0(e2,5^).(59) 

In the no-shear case and right-going waves (c = 1) one can recover (1250 : 
c = — 1 corresponds to left-going waves. The horizontal velocity to this order 
is 

Again, in order to invert (f59]) we have to specify u at a specific depth, 
(0 < zq < 1): n = n((^,r, zo). The result is 

(A 3 + 

r? = cu + e^-;^ -u^ — e(5^c^A'UA(', where A = — ; ^ -. (61) 

' 2(c2 + l) 6c(l + c2) 2c ^ ^ 

It is convenient to introduce a new dependent variable 
(which is NOT the Kodama transform) for which the equation is 

+ c2 + 1 



(c2 + l) 3c(c2 + l)^'^'^^ 



-e5 



■ . + c2 + 1 20*5 + 7c^ + 14c2 + 6 2a 
2cA 



C2 + 1 3(c2 + l)3 C(c2 + 1) 
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Next, we return to the original variables in ()63p . up to a scaling: T = 
^/et, X = ^/ex, see ([56]) : T = e~^T, X = ( + e~^CT, then 9^ = dx] 
edr = Ot + cdx ■ 

In conjunction, we add (ed'^ nVxxT — IJ'Vxxt)/{^ + c^) where fj. is an 
arbitrary real parameter. In the first term we substitute the leading order 
of 

(A A_ <i 

Vt ~ -Vx - £ , , , VVx 
+ 1 

and obtain 



Vt = -Vx 



c^ + c^ + 1 



VVx - eS"' 



1 — 3//C 
3c(c2 + 1) 



Vxxx + e<5^ 



C2 + 1 



[2cA(c2 + l) - 3^](c4+c2 + l) ^ 2c6 + 7c4 + 14c2+6 



2a 



(c2 + 1) 



^2 r2 
— e 



-c4 + 6c2 + 3 

3(c2 + 1)3 



3(c2 + 1)3 
+ c2 + 1 



VxVxx 



(c2 + 1) 



C(c2 + 1) 

VVxxx + 0{e\e6^) 



(64) 

The comparison between (I64p and ()36p gives the following possibilities 
for the parameters, cf. ([57|) - (f5^ : 



5(5 + l)(c'^ + 6c2 + 3) 

^ ~ 36(c2 + l)(c4 + c2 + l)' 

-26ci° + (3 - Ab)c^ + (21 - 66)0^ + (30 - 136)0^ + (27 - 26)c2 + 9 

66c(c2 + l)2(c4 + c2 + l)2 ' 

Therefore, according to the relation between zq and A in (|6ip . for a given 
propagation speed c, (136]) describes water waves at depth 



2 _ bc^^ + 86c^° + (186-3)c^ + (26b-21)c^ + (316-30)c^ + (126-27)c2 + 3b-9 

36(c2 + l)2(c4 + c2 + l)2 ■ 

In the case of right-moving wave without underlying flow (c = 1) we recover 
(pp . For the CH equation (6 = 2) this gives [cf. (Johnson 2003a)] 



_ V2ci2 _^ igglO _^ 33^8 _^ 31^6 _^ 32c4 - 3c2 - 3 

^/6(c2 + l)(c4 + c2 + l) ' ^ ^ 

and for the DP equation (b = 3) 

_ |c|Vcio + 8c8 + 17c6 + 19c4 + 21c2 + 3 

\/3(c2 + l)(c4 + c2 + l) ■ ^ ^ 
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The comparison between ([Ml) and ([l3]) gives [cf. (fiil) - (fi71) ] 



c'^ + 6c^ + 3 c2(c2 + 2)(2c'' + 6c2 + l) 

^ ~ 3(c2 + l)(c4 + c2 + l)' ~ "6c(c2 + l)2(c4 + c2 + l)2' 

Vc^2 _^ ggio _^ 18^8 _^ 26c6 + 31c4 + 12c2 + 3 

~ V^(c2 + l)(c4+c2 + l) ■ ^ ^ 

As expected, c = 1 in ([GT]) gives the result from (fiSI) . ~ 0.96. 

Without loss of generality we can assume that the underlying flow is 
propagating in the positive direction of the x-coordinate, i.e. A > 0. Then 
(|58p leads to the following restriction to the possible values of c: c > 1 (for 
the waves, moving in the direction of the flow, downstream) or — 1 < c < 
(for the waves, propagating upstream). The plot of the dependence of zq 
on c , according to the equations ([651) . ([66l) and (fHT]) is given on Fig|T]and 
Figl2j From Figl2]we notice that there is a region, where the function (j65p 
is not real: cq < c < 0, cq ~ —0.544, and therefore the CH (in this setting) 
is not a relevant model for these values of c. Although the graph of ()67p has 
a maximum zq ~ 1.0237 at c ~ —0.600. which violates the condition zq < 1, 
we notice that in the case of upstream propagation (Figl2]) for all possible 
values of c one can assume zq 1, i.e. equation ([1^ models very well the 
wave propagation on the surface {zq = 1) for all possible velocities. 



5 Conclusions 

In conclusion, CH, DP and (j42p describe in a direct way (without the use of 
the Kodama transform) the velocity (and, consequently all related variables) 
of shallow water waves at depths 0.71/i, O.SO/i and 0.96/i correspondingly (in 
the absence of a shear flow), where h is the depth of the undisturbed water. 
From a modeling point of view the advantage of the CH and DP equations 
over KdV consists also in the fact that they capture the wave-breaking phe- 
nomenon, cf. (Constantin 2000), (Constantin & Escher 1998) and (Zhou 
2004). These equations can also be used as water wave models in the pres- 
ence of an arbitrary shear flow. It is always more convenient to work with 
integrable equations, since their solutions are explicitly known or can be, in 
principle, explicitly constructed. For CH the solitons are stable patterns and 
thus physically recognizable, e.g. see the papers by Constantin & Strauss 
(2000, 2002), Constantin & Molinet (2001). The A^-soliton solution for CH 
is explicitly obtained by the Inverse Scattering Method in (Constantin et 
al. 2006) [see also the earlier works on the CH spectral problem by Con- 
stantin (1998, 2001); Constantin & McKean, (1999)]. In parametric form 
the A^-soliton solution for CH is obtained by: Johnson (20036), Li & Zhang 
(2004), Li (2005), Parker (2004, 2005a, 2005b), Matsuno (20056). Other 
types of explicitly known CH solutions are: multi-peakons (Beals et al., 
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Figure 1: Downstream propagation: plot of the dependence of zq on c. 

(f65]l for CH - solid line; ([Ml) for DP - dashed hne; ([67]) for equation (03]) 

- dash-dotted line. All cases have a horizontal asymptote zq ^ A= 0.577 

V 3 

as c — > oo. The CH graph has a maximum zq ~ 0.762 at c ~ 1.416; The DP 
graph has a maximum zq ~ 0.806 at c ~ 1.149. 




Figure 2: Upstream propagation: plot of the dependence of zq on c. ([65]) 
for CH - solid line; ((66]) for DP - dashed line; (f67]) for equation (113]) - dash- 
dotted line. The function ()65|) is not real for cq < c < 0, cq ~ —0.544; The 
graph of (|67|) has a maximum zq ~ 1.0237 at c ~ —0.600. 
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2003), periodic solutions (Gesztesy & Holden 2003), traveling- waves (Parkes 
&: Vakhnenko 2005) , (Lenells 2005a) . The constniction of multi-soliton and 
multi-positon solutions for the Associated Camassa-Holm equation using the 
Darboux/Backlund transform is presented in (Schiff 1998), (Hone 1999) and 
(Ivanov, 2005b). The A?^-soliton solution for the DP equation was recently 
derived by Matsuno (2005a), the multi-peakon solutions for DP arc obtained 
by Lundmark & Szmigielski (2003, 2005); the traveling waves - by Parkes &: 
Vakhnenko (2004) and Lenells (20056). There are similarities between the 
CH and DP equations, in a sense that they both are integrable and have a 
hydrodynamic derivation. However, it is interesting to notice that only CH 
has a geometric interpretation as a geodesic flow, cf. (Constantin &; Kolev 
2003), (Kolev 2004). 
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